HEXTE DETECTIONS OF HARD X-RAY TAILS IN SCO X-l 
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ABSTRACT 

We report the detection of a non-thermal hard X-ray component from Sco 
X-l based upon the analysis of 20-220 keV spectra obtained with the HEXTE 
experiment on-board the RXTE satellite. We find that the addition of a power- 
law component to a thermal bremsstrahlung model is required to achieve a good 
fit in 5 of 16 observations analyzed. Using PCA data we were able to track the 
movement of the source along the Z diagram, and we found that the presence 
of the hard X-ray tail is not confined to a specific Z position. However, we do 
observe an indication that the power law index hardens with increasing M, as 
indicated from the position on the Z diagram. We find that the derived non- 
thermal luminosities are ~ 10% of that derived for the brightest of the atoll 
sources. 

Subject headings: binaries:general — stars: individual (Scorpius X-l) — 
stars:neut r on — X-rays : st ar s 

1. Introduction 

Observations of hard X-ray emission from atoll sources (e.g. Barret, McClintock, & 
Grindlay 1996; Barret et al. 2000) have revealed the unexpected result that power law 
spectra extending above 30 keV may not be an exclusive signature of a black hole. Rossi X- 
ray Timing Explorer (RXTE) observations of hard X-ray emission from SAX J1808.4— 3658 
(Heindl & Smith 1998) demonstrated that hard X-ray flux is also present in other types 
of low-mass X-ray binaries (LMXBs). These observations refined the role of hard X-ray 
emission in distinguishing between neutron star and black hole binaries. 

Scorpius X-l is a high luminosity LMXB Z source, where the primary is a neutron 
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star with a low magnetic field. The presence of a non-thermal component in Sco X-l 
has been suggested (e.g., Peterson & Jacobson 1966; Riegler, Bolt, & Serlemitsos 1970; 
Agrawal et al. 1971; Haymes et al. 1972; Duldig et al. 1983) and the absence of such a 
component has also been reported (e.g., Lewin, Clark, & Smith 1967; Buselli et al. 1968; 
Jain et al. 1973; Johnson et al. 1980). In fact, high sensitivity searches have failed to detect 
such a component, placing strong upper limits on the non-thermal flux (e.g., Greenhill 
et al. 1979; Rothschild et al. 1980; Soong & Rothschild 1983). More recently, however, 
Strickman & Barret (2000) have reported the presence of a hard X-ray tail in Sco X-l 
using CGRO/OSSE. 

We used RXTE observations from the public database to search for hard X-ray tails 
in Sco X-l spectra. Sco X-l has been observed with RXTE several times, and our data 
sample (1997 April to 1999 July) contains 185,672 s of Proportional Counter Array (PCA) 
on-source time and more than 100 ks of on-source live-time in the High Energy X-ray 
Timing Experiment (HEXTE). We describe data selection and analysis in §2, results in §3, 
discuss the detection of a non-thermal component in § 4, and present our conclusions in § 5. 

2. Data Selection and Analysis 

In this work we used data from HEXTE to search for hard X-ray tails in Sco X-l in 
the ~ 20-220 keV interval, and data from the PCA to determine the Z position of Sco X-l 
in the color-color diagram. The PCA (Jahoda et al. 1996) is a set of 5 identical Xenon 
proportional counters covering the energy range 2-60 keV with a total area of ~ 7000 cm 2 . 
The HEXTE (Rothschild et al. 1998) consists of two clusters of 4 NaI(Tl)/CsI(Na) phoswich 
scintillation counters totaling ~ 1600 cm 2 of effective area in the 15-250 keV energy range . 
Each cluster rocks between source and two background fields to measure the instrumental 
background in near real time. The PCA background is estimated from particle rates and 
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the known diffuse X-ray flux. The PCA and HEXTE share a common 1° full width at 
half maximum field of view. Pertinent to this work, the HEXTE continuum sensitivity at 
100 keV, for ~5kseconds of live-time, is ~13mCrab (3a, scaled to in-flight backgrounds 
from the predictions of Gruber et al. 1996). 

We selected those observations from the RXTE public database (as of 2000 April) 
in which > 5000 s of HEXTE total on-source time is available, in order to achieve good 
sensitivity at high energies. Since some observations did not point directly at Sco X-l (in 
order to reduce the very high counting rate), we have also discarded observations for which 
RXTE pointing offsets were larger than 0.01°. These criteria resulted in selection of the 16 
observations displayed in Table 1. The resulting accumulated on-source HEXTE live-time 
is 104,238 s. 



2.1. Data Analysis 

We performed our analysis using HEXTE Archive Histogram Mode data, which are 
accumulated every 16 s, have spectral information in 64 compressed channels, and are 
present for all observations, independent of selected Science Mode. 

The data analysis was performed using software from FT00LS 5.0. For each 
observation, we divided the data into subsets according to position on the Z diagram 
(see below and Table 1). We then accumulated individual source and background count 
histograms for HEXTE clusters A and B and added them to form a single HEXTE source 
and background histogram for each subset. Response matrices for HEXTE A and B were 
generated and added (with weights of 0.57 and 0.43 respectively, according to their effective 
areas) . 

In order to be confident that any high energy excesses were real, we tailored our data 
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selection to optimize background subtraction. We verified that any remaining systematic 
errors as a result of our procedure were small and did not affect the results. We followed 
the standard practices extensively studied for HEXTE background subtraction by us and 
other workers (Rothschild et al. f998; MacDonald f999). First, we selected only data 
taken more than 15 minutes after satellite passage through the South Atlantic Anomaly 
(SAA). This allows any short-lived radioactive species produced in the detector material 
to decay and ensures that the systematic error of our linearly interpolated background 
estimate is <0.02% of the background. We verified this process by comparing - for each 
cluster - the spectra from the two separate background fields which are sampled alternately 
during observations. In each case, subtracting these spectra showed no excess emission, 
giving us confidence in the on-source background estimate. This also verified that there 
were no confusing sources in the background fields. We note, too, that none of the Sco 
X-l spectra show evidence for the prominent instrumental background lines at ~ 30, 70, 
160 and 190 keV (Rothschild et al. 1998) which we would expect to be apparent if there 
were over or under-subtracted backgrounds. Even though there was no obvious sign of 
residual background features, we allowed the normalization of the background spectra 
to vary in our fitting procedure. As expected, the fitted adjustment to the background 
level was always < 1 %. Even with this added degree of freedom, the significance of the 
detected high energy tails remained essentially unchanged. This also verified that live-times 
were calculated accurately, as systematic errors in the exposure would lead to incorrect 
background subtraction. Finally, at 100 keV (see Figure 1), the Sco X-l hard flux is of the 
order of 0.02 counts/s/keV, while the HEXTE background was typically 0.23 counts/s/keV. 
Thus, our measurements are at level of ~10% of the background. MacDonald (1999) 
showed that our background subtraction method, even for very long integrations, is free of 
systematic errors at a level of ~0.05% of the background. A more detailed discussion of 
HEXTE background systematics is given in Gruber et al. (1996). 
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We used the PCA data in order to constrain the position of Sco X-l in a color-color 
diagram. High resolution (~ 4 ms) lightcurves provided power spectra to characterize any 
quasi-periodic oscillations and/or continuum noise. Low resolution (16 s) lightcurves were 
used to generate color-color diagrams. We defined soft color as the ratio of the counts 
in the (5-7)/(3-5) keV energy bands and hard color as the ratio (7-20)/(5-7) keV. For 
observations 5/6, 12-15, 17/18 and 20/21 (see Table 1) the source was moving between the 
Normal and Flaring Branches. For those 5 dual-branch observations, we divided each into 
Normal and Flaring Branch subsets. This then resulted in 21 subsets. 

2.2. Spectral Fitting 

We used XSPEC 11.0 to fit the 20-220 keV HEXTE spectra with a combination of a 
thermal component plus a hard component. To fit the thermal component we tested the 
following models: thermal Comptonization (e.g. Lamb & Sanford 1979), blackbody, and 
thermal bremsstrahlung. The thermal bremsstrahlung model returned, in all of our subsets, 
the smallest x 2 > and it was utilized to represent the lower energy HEXTE flux. 

In the 5 subsets in which we detected the presence of a hard X-ray tail (see next 
section), a simple power law model, as well as ones with a spectral break, were evaluated. 
We found that the simple power law model described the hard component well, and it 
was used to represent the higher energy HEXTE flux. In all 5 subsets with significant 
non-thermal flux, a high temperature thermal component did not provide an acceptable fit. 

3. Results 

We developed two criteria to determine the presence of a hard X-ray tail in a particular 
spectrum: 1) a signal to noise ratio (SNR) > 5 in the 75-220 keV energy range, and 2) an 
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F-Test null significance for the addition of the hard component at a level of 10" 7 or less. 
We claim that we observed a hard X-ray tail only when both these criteria were fulfilled. 
We thus have 5 strong detections of non-thermal flux in Sco X-l spectra, with the hard 
X-ray tail extending at least to 220 keV. The spectral parameters derived for those subsets 
are shown in Table 2. Similarly we defined a strong non-detection when 1) the SNR in 
the 75-220 keV is < 1 and 2) the F-Test is at a level of > 10~ 5 . We thus have 4 strong 
non-detections. We note that several subsets formally have very significant F-Test values, 
but SNR values < 5. These typically result from deviations in the 30-50 keV range, possibly 
due to imperfect modeling of the thermal component. We cannot conclude that these 
subsets have truly non-thermal emission, although this possibility is open. In Figure 1 we 
show spectra both from one of our detections and one of our non-detections. 

The derived photon indices (r) for the 5 subsets where the tail was observed spanned 
a wide range: from —0.7 to 2.4, which differs from the OSSE results which have an average 
of 2.5 (Strickman & Barret 2000). It is suggestive that the two hardest photon indices were 
obtained when the source was on the flaring branch, and the softest were on the horizontal 
branch, showing a possible correlation between the hardness of the spectrum and the mass 
accretion rate M, which increases in the sequence HB ^NB -^FB (e.g., van der Klis 1996). 
We also note that both times the source was on the HB a hard X-ray tail was detected. 

The non-thermal 20-200 keV flux from the 5 detections varied by a factor of ~2.5 
from 0.63—1.56 x lO^ergscm" 2 s" 1 (see Table 2). Upper limits to the 20-200 keV flux 
from the 4 strong non-detections were estimated using a photon index of 1, the mean from 
the 5 detections. The one non-detection subset with equivalent sensitivity to the detection 
subsets has a 2a upper limit to the 20-200 keV flux of 1.0 x 10~ 10 ergscm~ 2 s -1 . This 
implies at least an order of magnitude in overall variability in the non-thermal component 
in Sco X-l. 
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In order to estimate the 2-20 keV luminosity, Ltfl^oi we use d a complex multicomponent 
model to heuristically fit the PCA spectra. We estimate the uncertainty in the derived flux 
due to uncertainties in the proper spectral shape to be ~2%. Using 2.8 kpc as the distance 
to Sco X-l (Bradshaw, Fomalont, & Geldzahler 1999), ranges over 2.1-3.2 x 10 38 

ergs s~ x for the 5 detection subsets. Similarly, the range of the 20-200 keV luminosities of 
the non-thermal component measurements is £20-200 = 5.9-15.0 x 10 35 ergss _1 . Thus the 
non-thermal component represents about 1 % of the total flux from Sco X-l. 

No correlation between the presence of the hard tail and the position of the source 
in the color-color diagram is indicated, i.e., we have observed the hard X-ray tail in all 
three branches. In addition the HEXTE thermal component is not correlated with the 
presence of the hard X-ray tail, i.e., the presence of the hard tail is indicated by neither the 
temperature nor flux of the HEXTE thermal component. 



4. Discussion 

From the lack of correlation between the presence of the hard tail and the position 
of the source in the color-color diagram source it would appear that the hard X-ray tail 
emission region is not associated with that part of the accretion disk which is believed to 
be responsible for the QPO behavior (van der Klis et al. 1996). These results also do not 
support a correlation between the presence of the hard tail and clustering in the color-color 
diagram, such as mentioned by Strickman & Barret (2000). We also note that this behavior 
is different from GX 17 + 2, in which a hard tail is correlated with the position in the Z (Di 
Salvo et al. 2000). Frontera et al. (1998) have also reported the detection of a hard tail in 
another Z source (Cyg X-2). 

Nevertheless we must point out that the observed photon indices may be related to 
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the position in the Z diagram. The hardest observed photon indices are in the FB and 
the softest ones in the HB. From this point of view the production of hard X-ray photons 
may be correlated with the accretion rate at the inner part of the disk. Motivated by 
the observation that Sco X-l is a variable radio source (e.g., Fender & Hendry 2000), 
we considered whether or not the hard X-ray emission could be explained in terms of 
synchrotron emission. While the measured photon indices are consistent with this idea for 
the HB and NB, the ones measured when the source was in the FB (consistent with 0) are 
not, since it would require an inverted distribution of relativistic electrons. Similarly, the 
HB and NB spectra could result from Comptonization, but the FB power law indices are 
again too hard. 

A difference between atoll sources and Sco X-l is also seen in the soft X-ray 
luminosity. The atoll source hard tails appear when the 2-20 keV luminosity is relatively 
low, < 10 37 ergss _1 ; whereas the soft flux in Sco X-l is essentially at the Eddington limit 
(~ 10 38 ergs s _1 ). It is also interesting to note that some models used to explain the hard 
X-ray emission of atoll LMXBs (e.g. Mitsuda et al. 1989; Hanawa, 1990) are applicable in a 
low luminosity state up to ~ 20 keV. This, however, is not the case of our observations. On 
the other hand, Sco X-l may be similar to some atoll sources in the luminosity of the power 
law component. The Sco X-l power law luminosity in the 20-200 keV range is comparable 
to the weaker atoll source hard X-ray luminosities reported by Barret, McClintock, & 
Grindlay (1996). 

Distinguishing neutron stars from black holes remains a key goal in high energy 
astrophysics. Recently it has been argued (Barret et al. 2000) that the luminosities in the 
1-20 keV (defined as L x ) and 20-200 keV (defined as L hx ) bands could be used to distinguish 
between black hole and low state neutron star systems displaying hard X-ray emission. 
In their argument, the critical luminosity, L crit = 1.5 x 10 37 ergss _1 , is the dividing line 
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between neutron star and black hole binaries with the neutron star systems emitting ~ 
L x when L x < L crit . Since Sco X-l is emitting Eddington-level soft fluxes, we cannot test 
that idea. On the other hand, our observed L hx does not contradict the idea pointed out 
by Barret et al. (2000) that only black hole binaries can have both L x and L hx above L crii , 
since our measured Lh x is below the L crit value (see Table 1). 



5. Conclusion 

We have observed significant detections of a non-thermal component in Sco X-l on 5 
out of 21 occasions with HEXTE on-board RXTE. We also have found that the presence of 
the hard X-ray tail is not associated with a particular position in the color-color diagram, 
but that the photon indices, instead, may be correlated with the position on the Z diagram, 
albeit with a very limited number of examples. We inferred that the non-thermal power 
law flux of the source varied by over an order of magnitude comparing our detections with 
the non-detections, and by a factor of ~ 3 in the 5 cases where the presence of the hard 
X-ray tail was strongly indicated. Our observed 20-200 keV luminosities are in agreement 
with only black hole binaries emitting hard X-rays with luminosity >1.5 x 10 37 ergs s _1 . 

This research has made use of data obtained through the HEASARC, provided by 
the NASA/GSFC. FD gratefully acknowledges FAPESP/Brazil for financial support under 
grant 99/02352-2. FD acknowledges Stefan Dieters, Wayne Coburn and John Tomsick for 
helpful discussions. We also thanks an anonymous referee for suggestions. This research 
was supported by NASA contract NAS5-30720. 
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Table 1: Selected RXTE data observations of Sco X-l 



OBSID 


MJD 


Nu. 


Z Pos. 


Live-time a 


SNR b 


F*e S i C 


Flux d 


Hard Tail (Y/N) 


20053-01-01-00 


50556 


1 


HB 


6341 


10.6 


4.4E-15 




Y 


20053-01-01-02 


50558 


2 


NB 


5434 


4.2 


1.7E-7 


4.05±8;|? 




20053-01-01-03 


50559 


3 


NB 


5896 


4.8 


3.9E-9 


„ . i n 75 
°- zu -0.74 




20053-01-01-05 


50561 


4 


NB 


1908 


0.8 


6.2E-5 


_ 11 77 

6.29±i.f? 


N 


20053-01-01-05 


50561 


5 


FB 


4593 


2.4 


7.7E-4 


^ . 1 QzL 

o 47+J--9 4 
-1.51 




20053-01-01-06 


50562 


6 


NB 


8558 


7.3 


3.5E-13 


_i_n qi 

8.63±8;|S 


Y 


10061-01-03-00 


50815 


7 


NB 


5484 


4.2 


2.7E-9 


^ 48 + 




20053-01-02-00 


50816 


8 


FB 


8145 


18.6 


5.8E-12 


13.6+J; , 


Y 


20053-01-02-03 


50819 


9 


NB 


5686 


3.5 


2.7E-9 


o.yi_ 95 




30036-01-01-000 


50820 


10 


FB 


6858 


5.8 


3.8E-7 


7-73±i;|l 


Y 


20053-01-02-04 


50820 


11 


FB 


4193 


4.9 


3.3E-4 


^ , .4.1 O^i 

°-^-1.00 




30036-01-02-000 


50821 


12 


NB 


4370 


2.2 


3.6E-4 


^•ou_ 90 




30036-01-02-000 


50821 


13 


FB 


2527 


1.8 


3.6E-4 


z - oo -1.06 




30406-01-02-00 


50872 


14 


NB 


2326 


2.5 


1.0E-6 


4-03±J;^ 




Qfi/in^ ni no nn 

oU4UD-Ul-Uz-UU 


OUo ( z 


1 K 


r r> 


olzo 


< u 


U.Uo 


^ o on 
< z.zu 


1ST 
In 


30035-01-01-00 


50963 


16 


HB 


6253 


6.2 


1.2E-14 


q fi q+0.92 

o.uo_ 88 


Y 


30035-01-03-00 


50965 


17 


NB 


3750 


1.3 


2.4E-8 


o on+0.86 
°- zu -0.89 




30035-01-03-00 


50965 


18 


FB 


2318 


< 


0.03 


< 2.67 


N 


30035-01-04-00 


50966 


19 


FB 


6230 


< 


0.38 


< 1.00 


N 


40706-01-01-000 


51433 


20 


NB 


8076 


2.8 


1.4E-5 


2 nq +0 - 55 

z.u»_ 53 




40706-01-01-000 


51433 


21 


FB 


2167 


2.2 


1.2E-3 


2 45+ 1 - 04 





Note. — HB=Horizontal Branch, NB=Normal Branch, FB=Flaring Branch; Nu.= observation number 
"total live-time for HEXTE, in s 
fc signal to noise ratio in the 75-220 keV energy range 
probability that a better fit occurred due to a random chance 

d flux in the 75-220 keV energy range, in units of 10~ 10 ergscm _2 s _1 ; apart from the 5 detections of a hard 
X-ray tail, the photon index used for the non-thermal component was our average value of 1 ; uncertainties 
are given at 90 % confidence level 
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Table 2: Observations of hard X-ray tail detections in Sco X-l 





Bremsstrahlung 


Power- 


-Law 




OBSID 


kT (keV) 


Flux a 


r 


Flux b 


X 2 (dof) 


20053-01-01-00 
20053-01-01-06 
20053-01-02-00 
30036-01-01-000 
30035-01-01-00 


4.8318;" 
4.5018:8^ 
4.3618:8^ 

a OO+0.03 
4 -^^-0.04 

A r. 1+0-08 
^• Oi -0.07 


q nq+o.39 

cr oq+0.30 
O.OO_ 27 

7.i6l8:iS 

9.63181? 

7.4518:8? 


1 7 c-+0.22 
l - ' °-0.20 

1 64+ - 29 

i -°^-0.27 

-0 17+0- 30 

u - i ' -0.33 

71 + - 63 

-U. i 1_ . 70 

9 07+O.33 
z -° ' -0.28 


i.56i8:i? 
o.9H8:io 

1 99+0.13 
1 - zz -0.14 

o.63l8:ii 
i.04l8:8§ 


36.76 (23) 
24.52 (23) 
29.15 (23) 
28.57 (22) 
24.01 (22) 



Note. — uncertainties are given at 90 % confidence level for the derived parameters of the model applied 
a Flux in 20-50 keV range in units of 10~ 9 ergscm^ 2 s _1 
fc Flux in 20-200 keV interval in units of 10 -9 ergscm~ 2 s _1 
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Fig. 1. — Spectrum resulting from a bremsstrahlung + power law fit to data subset 20053- 
01-02-00 (a), and the result of a bremsstrahlung fit to subset 30035-01-04-00 (b), showing the 
presence/absence of a hard X-ray tail in Sco X-l. Residuals are given in units of standard 
deviations (lower panels). In (b) the upper limits are 2a, including the 60-80 keV bin which 
experienced a — 3cr residual. The xl are 1-27 and 1.62 respectively. 



